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Abstract 

We reconsider the problem of even- visiting random walks in one dimension. This 
problem is mapped onto a non-Hermitian Anderson model with binary disorder. We 
develop very efficient numerical tools to enumerate and characterize even-visiting 
walks. The number of closed walks is obtained as an exact integer up to 1828 
steps, i.e., some 10^^^ walks. On the analytical side, the concepts and techniques of 
one-dimensional disordered systems allow to obtain explicit asymptotic estimates 
for the number of closed walks of steps up to an absolute prefactor of order 
unity, which is determined numerically. All the cumulants of the maximum height 
reached by such walks are shown to grow as k^^^, with exactly known prefactors. 
These results illustrate the tight relationship between even-visiting walks, trapping 
models, and the Lifshitz tails of disordered electron or phonon spectra. 
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1 Introduction 



Random walk is one of the most ubiquitous concepts of statistical physics. The large- 
time behavior of random walks is rather insensitive to details such as the presence of an 
underlying lattice or of an elementary time step, so that a vast range of them belongs to 
the universality class of Brownian motion, their continuum limit being described by the 
diffusion equation. 

Recently, in the context of surface growth, Noh, Park, and den Nijs [Q have introduced 
a model of random walks with the constraint that every visited site should be visited an 
even number of times. Such walks have been named even- visiting walks @, 0]. It has 
been argued that the non-local constraint of being even-visiting changes the universality 
class of random walks, with the typical extension of a walk of n steps scaling as n^/^'^+'^^ 
in dimension d, instead of the usual ra^/^. 

Derrida's argument for this behavior, to be recalled just below, demonstrates a 
deep analogy between the even-visiting walk problem, the trapping problem ^] of a 
diffusive particle which is absorbed by traps located at random positions, and the Lifshitz 
tails of the density of states of electron or phonon spectra of disordered solids near 
band edges. 

Consider for definiteness a random walk of n steps on the hypercubic lattice in di- 
mension d, with unit lattice spacing. The probability that this walk is confined within a 
sphere with radius R -C n^/^ scales as 

where j is the first positive zero of the Bessel function J(d-2)/2, so that j^/R"^ is the lowest 
eigenvalue of the Laplace operator in the sphere, with Dirichlet boundary conditions. 
Then, if ra S> R'^, most sites in the sphere must have been visited many times. Hence, 
to a good approximation, the probability that any given site has been visited an even 
number of times is 1/2, independently of other sites. Thus, among all walks that remain 
within the sphere, the fraction of even-visiting ones is about 




f{R) ^ exp(^-mn 2 R' 



where = n^^"^ /T(d/2 + 1) is the volume of the unit ball. Consequently, for large n, 
the main contribution to the number of even- visiting walks comes from 'Lifshitz spheres', 
i.e., optimal spheres whose radius i?L maximizes the product p{R)f{R), hence 

Both inequalities needed to justify the argument are satisfied by this solution, in any di- 
mension d. Eq. (|1 . 1| ) gives an estimate for the typical maximum extent of an n-step even- 
visiting walk from its origin. The total number of such walks scales as {2d)"'p{Ri^) f (Ri^) , 
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hence 

X{n) ~ (2ci)"exp -lt±(filn2)2/('^+2) . (1.2) 

In this paper we reconsider the one-dimensional case fl], ^, |^ in detail. We map the 
even-visiting walk problem onto a model with quenched disorder, analogous to a non- 
Hermitian Anderson model. Many concepts and techniques of one- dimensional disordered 



systems are then available p, y, |T0|. We develop efficient numerical tools to enumerate 
even-visiting walks. The parameter in the generating function of the walks governs the 
strength of disorder. The disordered system comes with an invariant measure, whose 
support is bounded for small disorder, and becomes unbounded at some critical point, 
at which the generating function of even-visiting walks develops a singularity. We first 
compute the escape probability ( |4.9| ), which is similar to the integrated density of states of 
the Anderson model. This allows to take advantage of the existence of a simpler problem, 
the random-mirror problem, for which many quantities can be computed analytically, and 
then translated in terms of even-visiting walks. Numerical checks show that the analogy 
between both problems indeed works at a quantitative level. We thus obtain in particular 
the estimate ( |6.11| ) for the number of closed even- visiting walks, and eqs. ([7.3|), ( [7.4| ) for 
the cumulants of the maximum height reached by such a walk. 



2 Definitions and mapping to a disordered system 

An n-step walk on a graph is a sequence voVi ■ ■ -Vn of vertices of the graph such that 
VqVi, . . . , Vn-iVn are edges of the graph. An even- visiting walk is such that for each vertex 
f 7^ fo, f„ of the graph, the number of indices i such that v = Vi is even. Note that this 
is a non-local constraint on the walk, so that it is not unlikely that statistical properties 
of even-visiting random walks are rather different from those of standard random walks. 

In this paper we shall concentrate on one- dimensional graphs, namely the segments 
with vertices 0, 1, . . . , or the half-line with vertices 0, 1, 2, . . . In each case, the edges 
connect nearby integers. 

We shall count even- visiting closed walks starting and ending at the origin (vertex 0). 
This problem can be reformulated as a one-dimensional disordered system ||^, Con- 
sider walks on a half-line with vertices numbered 0,1,2,... Each walk comes with a 
weight, which is computed as follows: each traversed edge contributes a factor and 
each arrival at site i contributes a factor Wi for i = 0, 1, 2, . . . We consider for the time 
being that t and Wq, Wi, . . . are commuting indeterminates. 

Let F(t, wq, Wi, . . .) be the generating function for walks P starting and ending at the 
origin, each walk being counted with its weight. Formally 

F{t, Wo, Wi,...)=}_^Vt l(7o • • • , 

p 

where S{P) is the number of steps, and Vi{P) the number of arrivals to site i. The 
only walk with steps has weight 1. Any other walk can be decomposed uniquely as a 
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succession of elementary blocks of the following type: a step from site to site 1, giving 
weight \/twi, a closed walk on the half line made of the sites 1,2, . . ., contributing to 
F(t,Wi,W2, . . .), and a step from site 1 to site 0, giving weight ^/iwQ. So 



F{t,wo,wi, . . .) = 1 + ViwiF{t,wi,W2, ■ ■ ■)ViwQ 

+ (^ViwiF{t,Wi,W2, ■ . .)Vtwoj +■ 

A resummation of the geometric series gives 

F{t, Wo, wi, ...) = ^ -, (2.1) 

1 - tWoWiF[t,Wi,W2, . . .) 

and leads to a continued-fraction expansion 

F{t,wo,w„...)= ^ l^^^ ■ (2.2) 

^ _ tWiW2 



These formulae are the starting point of all further considerations. Note that, in 
the expansion of F(t,wo,wi, . . .), the even-visiting random walks correspond exactly to 
monomials such that Vi{P), V2(P), . . . are all even. Then Vo(P) as well, because Vo(P) + 
Vi{P) + • • • is the number of steps, which is even for any closed walk. So, if the weights 
Wi are turned into independent random variables with vanishing odd moments and all 
even moments equal to 1, the average of ^(^,^0,^1, . . .), denoted by F{t,wo,wi, . . .), 
is the generating function for even-visiting walks on the half-line. In the same way, for 
finite segments, Wq, . . . ,wn are chosen as above, and wjy+i = 0. 

Now, a random variable with vanishing odd moments and all even moments equal to 1 
is simply a random sign: it takes values ±1 with probability 1/2. The above formulation 
can therefore be further simplified: if Wq, Wi, . . . wn are independent random signs, then 
so are = WqWi, e^^i = W1W2, . . ., ei = wn-iW^. 

The generating series of even- visiting walks on the segment 0, . . . , reads therefore 

FN{t,ei, . . . jEn) = :ir^ . (2.3) 

^ t£N-l 
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1 -tel 



It is clear from this expression that averages over the quenched disorder represented by 
the £i's will be even functions of t. As each power of t counts for two steps, the length 
(number of steps) of closed even-visiting walks is a multiple of four. 
Eq. (|2.3| ) implies the recursion formula 

Mt, ^1, • • • , Sn) = ^- ; r. (2.4) 
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If we write = '^n/'^n+i, with \l/o = = 1, then we have 



^iv+i - ^jv + te^^TV-i = 0. (2.5) 

This three-term hnear recursion relation defines a non-Hermitian one- dimensional An- 
derson model with ofi^-diagonal binary disorder. In this context, t can be interpreted as 



the strength of the disorder. The variables Fn, which obey the recursion relation (|2.4|) , 
are the associated Riccati variables |]10|, |1T|. 

It is clear that the t-expansions of F^(t, ei, . . . ,6^) and Fjv+i(t, ^i, ■ ■ ■ , £n) coincide 
up to order N (a closed walk of 2N steps cannot visit vertices higher than N), so all 
formula have a well-defined N —>■ oo limit. 

The t-expansion of Fn starts with 1, so its power F^ is well-defined as a formal 
power series in t, for a an arbitrary complex number, or even an indeterminate. Surely, 
only when a = 1 does this quantity have a simple interpretation as a counting problem. 
Nevertheless, considering arbitrary a also has its own interest (see section 8). Anyway, 
we set 

fN{t,a) = F§{t), fit, a) = lim fNit,a). 



3 Weak-disorder expansion 

If we formally expand F^{t) from eq. (|2.4| ) in powers of ^at, we find 
so that after averaging (and with m = 2n) 
and for N oo 

/(,a)^l + g/(,2„).-£|^, (3.) 

Let A/fc (resp. AfN,k) denote the number of even-visiting closed walks of Ak steps 
starting at the origin on the half-line (resp. on the segment [0, N]), and let us introduce 
the i-expansions 

hit, «) = E /iV,2n(«)t'", fit, a) = E f2nia)t'\ (3.2) 

n>0 n>0 

Eq. ( |3.1D yields the following recursion relation: 

r(a + 2n) 



fN,2mici) — X! fN-l,2m-2ni'2n) 



n=l 



r(a)(2n)!' 
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with initial conditions /Ar,o('^) = 1 for N > 0. The above formula is recursive both in 
and m, and it holds for > 1, m > 1. All the above quantities are non-decreasing 
functions of A^. Letting A^ — oo, we obtain 

- r(a + 2n) 

/2m(«) = E f2m-2n{2n) (3.3) 

n=i ^ {ajyzn). 

for m > 1, with the initial condition fo{a) = 1. 

This formalism can be used to compute the number of even-visiting walks in a very 
efficient way. To do so, it is better to trade the variable a for integers, to compute 
recursively an array /2m (2^*) forO<£<m<A; — 1, using 

m 

/2m (2^) = E /2m~2n(2n) 
n=l 

and finally to obtain the number of closed even-visiting walks with n = 4k steps on the 
half-line, starting at 0, as 

k 

A4 = /2fc(l) = E/2fc-2^(2^)- (3-4) 

The same idea can be implemented in the case of walks on a finite segment. 

With a few hours of computation on a workstation, using the MACSYMA software, 
Afk has been evaluated for k up to 457, and AfN,k foYN<k< 125. Note that AfN,k = A/fe 
for N > k. It is therefore possible to perform accurate investigations of the number 
of walks and of related quantities, such as the distribution of the maximum height (see 
sections 6 to 8). 

It is obvious that J\fk < 16'^, the total number of walks of 4/c steps. It is easy to 
generalize this (poor) bound to see that /2fc(«) < 2°'16''. So, the radius of convergence 
of f{t,a) in t"^ is at least 1/16. In fact, for real positive a, it is exactly 1/16, as can 
be seen from the above equations as follows. Suppose i > 1. From f(t,2i) = 1 + 
Emyi fit, 2m)t2™(^2£+2m-i^^ ^-^^^ t^-expansion of (l - t^^{*\~/)) fit, 2i) - 1 = 

i.2m.(2t+2m-l^ 



X^m^o,^^ ( 2m jfitj'im) has only non-negative coefficients. Hence the same is true 

of the if: -expansion of f{t,2e) - (1 - t^^{^^2e^)) ■ Inserting this bound for £ = m on 
the right-hand side of eq. (|3.1| ) shows that, for a > 0, the t^-expansion of fit, a) — 1 — 
Em^ot^"' (1 - ^^"'(^Tm^))"^r(" + 2m)/(r(a)(2m)!) has only non-negative coefficients. 

Note that the radius of convergence (in powers of t^) of (l — ^^™(^™^^)) very close to 
1/16 for large m. This implies that fit, a) is singular at = 1/16, giving the announced 
value for the radius of convergence. We have not been able to refine the above bounds in 
a useful way, starting from the weak-disorder expansion. In the next sections we obtain 
much better estimates by a different method. 
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4 Invariant measure and escape probability 



In order to get a feeling for the kind of singularities that appear at = 1/16, we first 
evaluate a quantity apparently unrelated to the counting of even-visiting random walks, 
namely the probability for the random variable F{t] ei, 62, . . .) to be larger than 2: 

Pesc(t) = Prob(F > 2). 

This quantity, referred to as the escape probability, vanishes identically for < 1/16, 
while it is non-zero for > 1/16. It is quite analogous to the integrated density of states 
in the Anderson model and similar disordered spectra p|, ^, |10|, |1^ . 



For the time being, consider a fixed real positive t, and define the (complementary) 
distribution function 

Rn{x) = Prob(F/v > x). 



The recursion relation (|2.4|) implies the following relation between the distribution func- 
tions of F/v and Fm-i- 

2Rn{x) = i?jv_i (^^) - Rn-i (i^) + Rn-1 {-j) - Rn-1 (^) + 20(-x), 



where 6 is Heaviside's the step function: 

e(x) = I 



for X < 0, 

1 for X > 0. 



For large N, Rn{x) approaches a limiting distribution function R{x), which defines an 
invariant measure dR{x), and obeys the Dyson-Schmidt equation [|l^, |Tl| 



The fractional linear mappings involved in this equation are the reciprocals of those 
involved in eq. 



4.1 The case < t < 1/4 

In this situation, both mappings involved in eq. ( p.4|) [or in eq. ( [4.1|) ] are hyperbolic. Let 
M{t) = (1 - VI - 4t)/(2t) be the smallest fixed point of the mapping a; i— > 1/(1 — tx), 
and let m{t) = 1/(1 + tM{t)). So, M{t) = 1/(1 - tM{t)) and -1/t < < m{t) < 1 < 
M{t) < 2 < 1/t. The mappings x 1/(1 — tx) and x ^ 1/(1 -|- tx) are respectively 
increasing and decreasing (on their intervals of continuity). Hence, if x lies in the interval 
I it) =]m(t), M(t)[, then the same is true of 1/(1 — tx) and 1/(1 -|-te). From -Fo(t) = 1, an 
inductive argument shows that -F/v(t, ei, . . . , ^tv) takes its values in I{t) for all N . Hence 
the support of the invariant measure is contained in I{t). In particular, we have R{x) = 
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for X > M{t) and R{x) = 1 for x < Tn{t), so that eq. (^4.1|) can be further simphfied for 
X G I{t) to 

The graph of R{x) is a (decreasing) devil's staircase, with plateaus at the dyadic numbers 
R = m/2", with n = 1, 2, . . . and 1 < m (odd) < 2" — 1. The support of the invariant 
measure dR{x), i.e., the closure of the set of points x around which R{x) is not a constant, 
is a Cantor set of measure zero. The heights of the plateaus do not depend on t, but 
their size monotonically increases with t. When t —>■ 0^, R{x) tends to the step function 
G(l — x). The other, more interesting limiting case t = 1/4 is investigated in next 
subsection. 



4.2 The case i = = 1/4 

In this borderline case, the mapping x i— >■ 1/(1 — t^x) is parabolic: it has a degenerate 
fixed point Xc = 2, which coincides with the upper bound of the interval /(1/4) = [2/3, 2]. 

The behavior of the distribution function R{x) near this upper bound will be our first 
example of an exponentially small, Lifshitz-like singularity. For x close to 2, eq. (|4.2| ) 
simplifies to 

2R{x) = R • (4.3) 

The above equation is actually an identity for x greater than 6/7, the image of the 
lower bound 2/3 by the mapping x H-i> 1/(1 + tcx). We set 

y = 77^, X = 2 - -, (4.4) 
2 — X y 

so that the mapping x ^— > 4(a;— l)/x corresponds to y y — 1. The general solution of 
eq. (|473|) therefore reads 

R{x) = 2-yA{y), (4.5) 

where A{y) is a bounded periodic function of its argument, with unit period. Periodic 
amplitudes are quite frequent in the realm of one-dimensional disordered systems ||^, |10 



p!3| . The present situation is analogous to that of the invariant measure of the Anderson 
model right at a band edge []T2[. Figure 1 shows a plot of the periodic amplitude A{y), 



obtained from exactly iterating the Dyson-Schmidt equation ( |4.2| ) a large enough number 
of times, starting from the initial data Ro{x) = 9(1 — x). 



4.3 The case i> 1/4 

In this situation, the mapping x 1/(1 — tx) is elliptic: its two fixed points acquire an 
imaginary part. With the parametrization 

^ = (0 < </. < n/2), (4.6) 

4 cos^ (p 
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the fixed points read x± = 1 + exp(±2i</)). As a consequence, the support of the invariant 
measure is unbounded. A fraction of the measure escapes above 2, hence the name 
'escape probability' for the quantity Pgsc = -R(2). 

When t — 1/4 ?s (p'^/A is small, the fixed points x± ~ 2 ± 2i0 are close to Xc = 2, 
so that the escape phenomenon takes place through a narrow channel. It is therefore 
legitimate to deform the result ( [4. 51) as follows. Instead of eq. (|4.4|), we set 

e-^ = ^^^, x = 2cos0^^^^M^ (0<^<vr), (4.7) 
X — X- sin yj 

so that the elliptic mapping x (x — l)/{tx) corresponds to ^lJ ^-^ ip — (f). The angle 
ratio is therefore the appropriate deformation of the co-ordinate y (up to an additive 
constant). There is actually no such constant, as eqs. ( [4.4| ) and ( |4.7| ) yield ip = y(j)+0{(f)^). 
Eq. ( [4.5| ) therefore implies 

/?(x) ^ 2-^/'^A(V^/0), (4.8) 

in the regime where is small, while the angular variable ip is arbitrary. In particular, 
X = Xc = 2 corresponds to ip = 7r/2 + 0, hence the following prediction for the escape 
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probability: 



Pesc - 2-'^/(^^)- A(vr/(20)) ^ 1 exp -^^^ A 1 , (4.9) 

^ ^ ^ \ Vt-1/4; V4 - 1/4/ 



which is the resuh we were after. The escape probabihty is therefore exponentially 
small in {t — l/4)~^/^, with an asymptotically periodic amplitude A, already entering 
eq. ( [4. 51 ), and shown in Figure 1. The result ( [4. 91 ) is fully analogous to the Lifshitz 
tail of the integrated density of states in the Anderson model and similar disordered 



spectra M, 10, 12 



We have not been able to obtain by direct means estimates of the (more interesting) 
singularities that govern the asymptotics of the number of even-visiting random walks. 
However, the result ( [4.9|) for the escape probability is quite suggestive, especially because 
of its formal analogy with the Lifshitz tail for the integrated density of states of elec- 
tron or phonon spectra. The crucial ingredient in this result is that the transformation 
X (— > 1/(1 — tx) turns from hyperbolic to elliptic at t = tc = 1/4. Equivalently, the 
singularity ( |4.9|) is due to the occurrence in the recursive formula ( |2.4| ) of regions made 
of a large number, of order 7r/(20), of consecutive positive e's, which are one-dimensional 
analogues of Lifshitz spheres. The other transformation x ^ 1/(1 + tx) has no impact on 
the exponentially small form of the singularity. Its influence is only felt on the periodic 
modulation. 

In the next section, we introduce a simpler model, the random-mirror model, which 
keeps the essential feature of the original even- visiting walk problem, i.e., the transfor- 
mation a; I— >■ 1/(1 — tx), but for which quantities of interest are directly computable by 
elementary means. The random-mirror model is similar to the binary random harmonic 
chain where a finite fraction of the atoms has an infinite mass, first considered by Domb 



et al. |T^. In this limiting case, the system splits into an infinite collection of independent 
finite molecules, so that the integrated density of states, among many other quantities, 
can be evaluated by simple algebra. The same simplification occurs for diffusion in the 



presence of perfectly absorbing sites [T^. By analogy with these situations, we claim that 



the random-mirror model and the even- visiting walk problem have similar exponentially 
small singularities. The simple and explicitly computable periodic amplitudes of the first 
model are just replaced by (admittedly complicated) unknown periodic functions. In 
the leading terms of large-order t-expansions, this replacement is just responsible for an 
absolute prefactor. 



5 A case study: the random- mirror model 

If the distribution of the variables is modified, so that they take value or 1 with 
respective probabilities p or 1 — p (before they took value —1 or 1 with probability 1/2), 
the equations obtained previously for even- visiting walks are modified in a straightforward 
way. We assume p 7^ 0, 1. The physical interpretation of the present model is that any 
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site i = 1,2,... is totally reflecting with probability p, so that the particle freely moves 
between the origin and the first reflecting site, hence the name 'random-mirror model'. 

From a more mathematical viewpoint, this amounts to replacing the map x ^ 1/(1 + 
tx) by the constant map x ^ 1, without changing the second map x i— > 1/(1 — tx). We 
shall also consider the case of a general constant map x ^ X (instead of x t-^ 1), where 
X is an indeterminate (most of the time specialized to a real number X < 2). 

We shall compare in detail the random-mirror model and the even- visiting walk model, 
so that it is useful to change notation: for the random-mirror model we use G, g, and S, 
instead of F, f, and R. 



5.1 Weak-disorder expansion 

Let us first fix X = 1. From the recursion 

GN{t 



1 — ^^AfG'Ar_i(t, £1, . . . , ^TVi) ' 

we get 

T(a + n) 



gjvit, a) = G%{t) = 1+pY^ gN-i{t, n)t 



n>i r(a)n! 



which in turn leads to a recursion relation for the coefficients of the expansion gN{ti ot) 
Y.n9N,n{Q)t'', namely 

^ V{a + n) 



n=i r(Q;)n! 
with initial conditions gN,o{ci) = 1 for X > 0. For N oc, we get 

g{t, a) = l+pY: 9{t, r^)r ^^" + ^\ (5.1) 
^1 r(a)n! 

and 

- J{a + n) 

gm[a) =P2^ gm-n[n) ^ (5.2) 

for m > 1, with the initial condition go{a) = 1. We can repeat the argument given in 
the previous section for fit, a) to show that for positive a, the t-expansion of git, a) 
has radius of convergence 1/4. The above formalism is easily extended to generic values 
of X. For example, eq. ( p.l[ ) becomes 

git,a) = l-p + pX^+pY.g{t,n)t- ^ ^ ' 



„>i r(a)n! 



while eq. ( |5.2| ) remains unchanged, with the initial condition being modified to go{a) 
l-p + pX''. 
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5.2 Escape probability 



We now evaluate the escape probability by more direct means than for the case of even- 
visiting walks. The random variable G is equal to Qn with probability (1 — p)p^, where 
the Qn are defined recursively as follows: 

This definition gives a meaning to Qn for any integer n, positive of negative, and shows 
that Qn is a linear fractional function of X. The distribution function S{x) = Prob(G' > 
x) reads explicitly 

S{x) = {l-p)J2p^e{Qn-x). (5.4) 

n>0 

For the random-mirror model stricto sensu, we have X = 1. Suppose t > 1/4. With 
the parametrization ( |4.6| ), we infer from the recursion (O) the explicit expression 

^ / cosfn + 1)</)\ 

Q„ = 2 l-sin0 . n . 5.5 

y sm(n + 2)0 J 

Then Qn > 2 is equivalent to either sin(n-|-2)0 > and cos(n+l)<^ < 0, or sin(n+2)</) < 
and cos(r;, + 1)0 > 0. The smallest n such that Qn > 2 belongs to the first case, and 
readsQ uq = [7r/(20)] — 1. When t is close to 1/4, i.e., is small, many n's following 
no, up to Til ~ 2no, have Qn > 2. As a consequence, Pesc = >S'(2) can be estimated as 
(1 -p) Eno<n<nii'", leading to the result 

Pesc^P^-^^""^^-', (5.6) 

up to a correction of order 0{p'^^'^). 

When the initial point is not 1, but a generic X, eq. (p. 51) generalizes to 



/ . ^ cos(n + 1)0 + (1 - X) cos(n - 

Q„ = 2 l-sm0 — — 77^^ 7 — . 5.7 

Y sm(r;, + 2)0 + (1 — A ) sm n0 J 

The smallest n such that Qn > 2 reads no{X) = [7r/(20) — A(0, X)], with 

1 / X \ 

A(0,X) = -arctan I tan0) . (5.8) 

V2 — X / 

Hence 

Pesc(x)^prv(2^)-^(*-^)i, 

again up to a correction of order 0{p^^'^). For small 0, this result can be recast in a form 
similar to eq. ( [4.9| ), namely 

Pesc(X)^p-/(2*)Ax(7r/(20)), (5.9) 



^The symbol \x~\ denotes the smallest integer larger than or equal to x, while lx\ denotes the largest 
integer smaller than or equal to x. For x not an integer, we have \x~\ — \_x\ + 1, and the fractional part 
of X is defined as Frac(x) = x — [ccj = x — \x] + 1. 
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where the periodic amphtude Ax{x) has the exphcit expression 

Ax{x) = p-F'-ac(a;-Ao(X))+l-Ao(X)^ 

with 

Ao(X)=A(0,X) = ^. 

This resuh shows that the random-mirror problem at p = 1/2 and even- visiting 
walks have the same leading singularities (compare eq. (|5l9| ) and eq. ([4.9|) ): the escape 
probability of the random-mirror problem (X = 1), of the generalized one {X generic), 
and of even-visiting walks only differ by their periodic modulation. The same property is 
known in the case of spectra of disordered systems: details of the distribution of random 
masses, random site potentials, and so on, only enter the shape of the periodic amplitudes 
of Lifshitz tails g |ig, |12[. 



6 Number of even-visiting walks 



The purpose of this section is to establish the asymptotic estimate ( |6.11|) for the total 



number of even- visiting walks starting and ending at the origin, using eqs. (|3.2| ) and (|3] 

The escape probability Pesc was easy to estimate in the even-visiting walk problem, 
because it is a purely singular quantity, which vanishes as the parameter t approaches 
tc = 1/4. To the contrary, the singularities of the functions f{t,a) or g{t,a) at t = 1/4 
are less easily grasped, because the latter quantities also have a regular part. 

We shall first consider the random- mirror model with X = 1. In this case, we shall 
see that there is a simple relation between the discontinuity of g{t, 1) for t > 1/4 and the 
escape probability. A similar correspondence is expected to hold for the discontinuity 
of /(t, 1), up to a periodic modulation. These discontinuities provide a direct way to 
estimate the large-order behavior of t-expansions. The latter will lead us to estimate the 
number of even- visiting random walk (eq. ( |6 . 1 1| ) ) . 

6.1 Discontinuity of g{t,l) 

Let us consider first the random- mirror model in the simpler situation where X = 1. By 
analogy with eq. (^.4|), we have the closed formula 

g{t,a) = {l-p)Y.P''Qn- (6.1) 

n>0 

Their recursive definition ( |5.3| ) shows that the Qn{t) are rational functions of t. Their 
explicit expression (|5.5|) yields the decomposition 



4 L(-+i)/2j / ^ \ vrm 
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It is therefore clear that all the poles of the rational function Qn{t) lie on the half-line 
]l/4, +oo[. We normalize the discontinuity of a real-analytic function f{t) on the real axis 
as Df(t) = (=Fl/7r) Im/(t±iO), in such a way that the discontinuity of 1/t is exactly S(t). 
Let t be real and such that t > 1/4, and consider the integrated discontinuity 

D{n,t) = f DqM^u. 

Jl/A 

Eq. (§3) yields 

. L("+1)/2J / . X 

D(n,t) = — y Qlt tan^ 

^ ^ n + 2 \ 4cos2^y 



n + 2 



(6.3) 



Eq. ( |6.1| ) implies that the integrated discontinuity of g{t,l) reads 

D,{t) = {l-p)J2p-D{n,t). 

n>0 

Eq. ( |6.3|) leads to the following explicit result, in terms of 0: 



1 n L("-1)/2J . s 



2 vrm 
tan . 



which we reorganize as 



2 vrm 

^ ^ tan . 

P m>ln>2m-l \ ^ / ^ 

Finally, because < < 7r/2, the condition ixm/n < is always more stringent than 
n > 2m — 1, hence 

m>l n>\-Km/4i\ 

For small the contribution of m = 1 is exponentially larger than the other ones, so 
that the asymptotic behavior of -Di(0) is 



0' 



3 



^ -^pl^/'^l-^ (6.4) 
vr 

The above formula entirely comes from the pole of Qn{t) closest to the critical point 
tc = 1/4. This gives the clue for treating the case of a general constant map. In that 
case, Qn, as given by eq. (p.7|) , has poles for tan(?2 + 1)0/ tan0 = — X/(2 — X). Let 
(f)i{X,n) be the smallest solution of the latter equation. One checks that, for large n, 
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The contributions of the closest poles for each n to the total integrated discontinuity are 

_n - AV^ si^^^'^i n cos(n + 1)01 + (l-X)cos(n- 1)01 
^ ^^4^cos3 0/ (n + 2)cos(n + 2)0i + (l-X)ncosn0i 

This expression looks formidable, but it simplifies drastically at small 0. From the 
definition of 0i(X, n), the smallest n giving a non- vanishing contribution to the above 
sum is such that 

tan(n + 1)0 ^ X tann0 
tan0 ~ 2 — X tan0 ' 



so that n = [vr/0 — A(0, X)] — 1, with the definition (|5.8| ). Hence 

a3 



Dx{(j)) ~ p 

71 



I.e., 



with 



D^(^)^-^p-/^Bx{7c/(P), (6.5) 



Bx(x) = p-F-(..-Ao(X))-Ao(X) ^ (6.6) 

p 

The exponentially small factor p'"^'^ in the result (|6.5|) for the discontinuity of g{t, 1) is 
just the square of the corresponding factor in the result (|5.9| ) for the escape probability. 
This correspondence, which holds for any value of X, is expected to hold as well for the 
even- visiting walk problem. The simple relationship between the periodic amplitudes Ax 
and Bx is, however, a peculiarity of the random-mirror problem, that we do not expect 
to hold in general. We quote for further reference the constant Fourier component of the 
periodic function Bx{x), 

Bf = [' Bxix) dx = p-Ao(X)-l ^ }_ZP p-2/(2-X)_ 7) 

JO |lnp| |lnp| 



6.2 Asymptotic results 

We again consider first the generalized random- mirror model, with arbitrary X. The 
knowledge of the integrated discontinuity Dxif) allows to compute the asymptotics of 
the t-expansion of g{t, 1). By definition, the coefficient of order k in this expansion reads 

..(1) = i^9{tA) = = - r''(4cos^0)^^^^d0. (6.8) 

' J 2Tiit^+^^^ ' ' A/4 t^+^ dt Jo ^ d0 ^ ^ ' 

For large k, the last integral is dominated by small values of 0, where cos0 can be 
expanded to quadratic order and exponentiated, while the asymptotic result ( |6.5|) holds 
for Dx{4>)- It turns out that only the constant Fourier component B^^^ of the periodic 
amplitude Bx{x) matters (see below eq. ( |6.10D ). We thus obtain 

gJl) ^ \\np\ Bf 0e-^'^'-"l'"Pl/'^ d0. (6.9) 

Jo 
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This integral is then evaluated by the saddle-point approximation, with the saddle-point 
value of (j) being (jy^ = (vr |lnp| /(2k)Y^^. We are thus left with the explicit asymptotic 
result 

^,(1) ^ 4 bP (^)'^'exp (-^ {2n\lnpr kf'^ 4^ (6.10) 

which only depends on X through (see eq. (|6.7|)). 

It is now clear why only the constant component B^^ matters, just as in trapping 
problems p, |10|- Consider the Fourier expansion Bx{x) = Efc^j^^e^-i^^ The con- 



tribution of the Fourier component Bx^ is obtained by replacing in eq. (|6.1CI| ) |lnp'^'^^ 



ATk ~ 5(°)(21n2)^/3 /^y/%^p ((2vrln2)2 A;)'^') (6.11) 



by (— Inp — 27rifc)^/^, a number with a strictly greater real part, for any 7^ 0. The 
oscillations of the result ( |6.5|) therefore get damped by exponentials of fc^^^ in the large-/c 
asymptotics of quantities such as (7^(1). 

Coming back to the even-visiting walk problem, using the correspondence between the 
singularities at t = 1/4 of f{t, 1) and of g{t, 1) with p = 1/2, we arrive to the following 
asymptotic expression for the number A4 = /2A:(1) of closed even- visiting walks of 4k 
steps: 

J J ~^""V 2 

Note that the case of even- visiting walks has the symmetry t ^ —t, so that f{t, 1) 
has two cuts, at ] — cxD, —1/4] and [1/4, -|-cxd[, which cancel each other at odd orders in t, 
and add up constructively at even orders. 

In the expression (|6.11|) , B^^^ stands for the constant component of the periodic 
amplitude B{x) of the even- visiting walk problem, which cannot be predicted analytically. 
The exact enumeration procedure described in section 3 has been carried out up to 
k = 457, i.e., 4fc = 1828 steps. The number of these walks (exactly evaluated as an 
integer) reads ACtsr ~ 3.36829575 x 10^^^. Figure 2 shows a plot of the ratio A4 to its 
asymptotic behavior J\f^^\ defined as the result ( |6.11| ) without its prefactor B^^K The 



data exhibit a smooth convergence in This form of the leading finite-fc correction 

is indeed expected for all asymptotic estimates, as it corresponds to regular corrections 
to leading critical behavior, of relative order t — t^, via the expression of the saddle-point 
0c given below eq. ( |6.9| ). We thus obtain the accurate estimate 

^ 1.760. (6.12) 



7 Distribution of the maximum height 

The purpose of this section is to investigate the distribution of the maximum height 
reached by an even-visiting walk of 4k steps, starting and ending at the origin. 

The generating function for even- visiting walks with maximum height M is the differ- 
ence /m(^, 1) — fM-i(t, 1), and its coefficient of order 2k is /M,2fc(l) — /Af-i,2fc(l)- Hence 
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Figure 2: Plot of the numbers A4 of even-visiting 4A:-steps closed walks, divided by their 
asymptotic formM^^^ defined in the text, against up to A; = 457. Full line: least-square 

fit of data with k > 50, with intercept = 1.760. 

the probability that the maximum height in an even-visiting walk of 4k steps be M is 

/Af,2fe(l) — /Af-l,2fc(l) 



n 



M,k 



f2k{l] 



(7.i; 



The generating function of the cumulants of M reads therefore 



j>l 3- \M>0 



In 



1 -e" 



, ... E/M,2.(l)e^^^ . (7.2) 



Let us again make a detour through the random-mirror model. It is straightforward 
to extend eq. ( fr2|) to the latter situation, by simply replacing /'s by g's, and Z by Y, 
even if the probabilistic interpretation of the formula thus obtained is delicate. 

The identity 



TV 



gN{t,l) = {l-p)J2p'^Q^it) 

n=0 



implies 



Y,(z) = In 



Ck{p) 
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with 

The large-fc estimate ( |6.1U| ) for gk{l) leads to 



Coming back to the even-visiting walk problem, using the correspondence between 
both problems for p = 1/2, we obtain 

Zkiz) ^ ^ {An'ky/-' ((In 2)2/3 - (In 2 - z)^^^) . 

Hence all the cumulants of the distribution of the maximum height M of 4A;-step walks 
are asymptotically proportional, namely 

<M^>^a,fc^/3, (7.3) 

with explicit prefactors 

"^ = ^r(l/3f (2^)'^^a"2)2/3-^ (7.4) 

i.e., ai = 3.847495, as = 1.850254, = 3.559136, = 11.981080, and so on. 

As the number of steps becomes large, the distribution of the maximum height M 
therefore gets more and more peaked around its mean value ai fc^/^. The bulk of this 
distribution is a narrow Gaussian, with standard deviation a2^'^k^^^. This result is in 
sharp contrast with usual random walk, for which the largest positive extent M scales 
as t^/2, and the ratio ^ = M/t^^"^ admits a non-trivial limiting probability law |15 . 



The above analytical results have been checked against data obtained from the exact 
enumeration procedure of section 3. Figure 3 shows a plot of the first four reduced 
cumulants -C ^ (j = l;---;4). The clear linear behavior in k^^^, with unit 
slope, demonstrates a quantitative agreement between the exact data and the asymptotic 
results (pD, (|73). 



8 Amplitude ratios 

Up to now, we have only considered averages of Riccati variables with respect to the 
quenched disorder, such as f(t, 1) and g{t, 1), which have a direct combinatorial interpre- 
tation. However, the even-visiting walk model and the random-mirror model also display 
an interesting dependence in a, which will be investigated in this last section. 

To start with, let us consider the random-mirror model with X = 1, and assume that 
a = £ is a positive integer. We can use eq. (|6l^ ) to obtain the pole structure of Qn(t), 
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5 




Figure 3: Plot of the first four reduced cumulants <C S> /aj of the maximum height M 
of 4A;-step even-visiting walks, against k}/^, up to A; = 125. The full lines with unit slope 
demonstrate the quantitative agreement with the analytical results (|7.3| ), (^^). 



which governs the asymptotics of its t-expansion. As before, it is clear that the pole 
closest to the critical point (m = 1) gives the leading contribution. So, we write 



Q'nit)=Y. , '' ' ' (8.1) 

1 -4tcos2^) 



n+2 J 

where the dots stand for a part which is regular at t = 1/(4 cos^^^). It turns out 
that the coefficients Cj{n,l) have a simple behavior when n is large. Indeed, let us write 
(j) = 7r/n-7rX/n^ + 0{l/n^), so that t = l/(4cos2 ^)-7r2 A/(2n^) + (9(l/n^), and re-express 
Qn-2{t) in terms of A. Using eq. ( |5.5| ), we can check that Qn-2 = 2(1 + 1/A) + 0(l/n), 
so that 



On the other hand, a direct substitution of the above expression for t into eq. ( |8.1| ) yields, 
for large n, 



• n \ W 



19 



Comparing the above two formulae, we obtain, again for large n, 



Cj{n,i) 




Now, the t-expansion of eq. (|8.1|) is straightforward. A pole of order j gives an extra 
factor {j + k — l)!/((j — 1)!A;!) compared to a pole of order 1. The sum over n is then 
performed by the saddle-point approximation, with the saddle-point value of n being 
ric ^ (27r^A;/ |lnp|)^/^. We are thus left with the following explicit values of the amplitude 
ratios: 

^' »-»9t(i) p,(3-iy\j) ^ ' 

It turns out that this result is independent of X, but we shall not go into any more detail. 

Quite unexpectedly, the amplitude ratios are also determined as the solution of 
an eigenvalue problem. The starting point is equation (|5.2|) , which we divide by gm{^)- 
Because the radius of convergence of the t-expansion of g(t, a) is 1/4, we expect that, for 
large m and fixed n, gm-n{n) ~ 4~"(7m(n). So, we get 



n>i 

In particular, fixing a = £, a. positive integer, 

e(,).pge(„)4-"(^i^. (8.3) 

This equation implies that the (infinite) matrix A4 with positive entries Ai^^n = + 
n — l)!/(4"(£ — l)!n!) for i,n > 1 has an eigenvector with eigenvalue 1/p and positive 
components Q{i), for every p g]0, 1]. We have therefore an explicit example of why 
the Perron-Frobenius theorem (stating that a finite positive matrix has exactly one 
eigenvector with positive components) cannot extend to infinite matrices without further 
hypotheses. Note that the matrix A4 is highly asymmetric, has unbounded entries, and 
a divergent trace. 

It is interesting to compare the linear system (|8.3| ) with eq. (|5.1|) at t = 1/4, which 
reads 

^^(1/4, «) = 1 + p E ^^(1/4, 

^1 r(«)n! 

In particular, fixing a = i, a. non-negative integer. 



(7(1/4, i) = (7(1/4, 0) + p E (7(1/4, n)4- 



n>l 



{e-l)\n\ 



So, despite the fact that analyticity for \t\ < 1/4 and continuity at t = 1/4 fix g{l/4:,a) 
unambiguously, the are a (positive) ambiguity of the solution of eq. (|5.1j ) restricted 
to t = 1/4. This is again incompatible with the Perron-Frobenius property. 
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Because of these strange properties, it is reassuring to show that the exphcit for- 
mula ( p.2| ) for the Q{i) indeed satisfies (^]3|), which was derived only heuristically. The 
proof relies on the following identity: 

y Si±3Ls^ = e^£' y 

which holds for i a positive integer, and can be checked by an explicit expansion of its 
right-hand side in powers of s. We thus obtain the explicit formula 

When a is large and positive, at fixed p 1, the sum over j can be evaluated by the 
saddle-point approximation, the saddle-point value of j being ~ (|lnp| a)^/^ + |lnp| /2, 
yielding 

^(«) ^ \ {-J^) «)'^') 2"- (8-5) 

In the present case of the random-mirror model, we have been lucky enough to obtain 
the exact amplitude ratios (|8.2| ). Had it not been the case, eq. (|8.3|) would provide a 
convenient tool to get accurate numerical values for these ratios. In fact, this alterna- 
tive approach leads to a better understanding of how the Perron-Frobenius theorem is 
bypassed. If the matrix J^e^n is truncated to n < £max, then, whatever £max is, the 
spectrum does not contains any eigenvalue larger than unity. So, the right trick is to 
truncate Aii^n to £ < imax, n < ^max + 1, to fix the normalization = 1, and to solve 
for 6;(2),...,'6;(£^ax + l). 

Coming back to the even-visiting walk problem, we shall follow the latter strategy. 
We expect, by analogy with the previous situation, that = limfc^oo (/2/c(«)//2fc(l)) 
are well-behaved limits. Then, starting from eq. ( p.3| ), we can repeat the above argument, 
to show that 

^1 r(a)(2n)! 
and, fixing again a = i, a non-negative integer, 

^(«)^g^(2„)16-"|±^^. (8.6) 

This result implies that eq. (|3?T|) that determines fit, a) is ambiguous at = 1/16. By 
truncating the system ( ^.6|) as explained above, we can obtain the amplitude ratios J-'{i) 
with very high accuracy, namely J^(2) = 5.591806, J^(3) = 22.285850, J^(4) = 77.059126, 
and so on. 

These predictions have been checked against exact data for the ratios /2fc(^)//2fc(l), 
obtained by means of the enumeration approach of section 3. Figure 4 demonstrates a 
quantitative agreement, for £ = 2 to 4. The data smoothly converge to the predicted 
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limiting values J-'{i), with (small) A;~^/^ corrections. This reinforces our assertion that 
even-visiting walks and random-mirror models have very similar asymptotic behaviors, 
up to multiplicative prefactors. 
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Figure 4: Logarithmic plot of the ratios /2fc(^)//2fc(l) against /c^^/^, for £ = 2 to 4, and k up 
to 400. Full lines: least-square fits of data with k > 50. The intercepts match very accurately 
the predicted limiting amplitude ratios (symbols on vertical axis). 



9 Discussion 

In this paper we have reconsidered the problem of even-visiting random walks, and ob- 
tained several kinds of exact or asymptotic results in one dimension. 

We have mapped the even-visiting walk problem onto a non-Hermitian Anderson 
model ( |2.5|) . The weak-disorder expansion of the latter model provides very efficient 
numerical tools to enumerate and characterize even-visiting walks. We have thus been 
able to evaluate exactly, among other quantities of interest, the total number A4 of closed 
even- visiting walks up to = 457, i.e., 4fc = 1828 steps, going thus far beyond previous 
works. Indeed, our result is to be compared with the exact result 0] up to 80 steps, and 
with the approximate numerical simulation |l|] up to around 1000 steps. 

The mapping to the non-Hermitian Anderson model (|2.5| ) also makes many concepts 
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and techniques of one- dimensional disordered systems available 10|. The analogy with 



Lifshitz tails, which was already clearly apparent from Derrida's argument ^ in any 
dimension, has been corroborated at the level of analytical tools in the one-dimensional 
situation. The escape probability, investigated in section 4, is fully analogous to the 
integrated density of states in disordered spectra ^ [1^, |T2[, with an exponential 



behavior, modulated by an oscillatory amplitude, which is asymptotically periodic in 
the relevant variable. Periodic amplitudes are actually ubiquitous in one- dimensional 



disordered systems, when the quenched disorder has a discrete distribution g, |T0|, |T3 . 

In order to investigate more complex quantities, such as the generating functions of 
even-visiting walks, we make an extensive use of the random-mirror model. The latter 
model is similar to the random harmonic chain where a finite fraction of the atoms 
has an infinite mass, first considered by Domb et al. [Q: quantities of interest are 
directly computable by elementary means. The results thus obtained still hold true for 
the even-visiting walk problem, the only difference being that the simple and explicit 
periodic amplitudes of the first model are replaced by unknown periodic functions. The 
same phenomenon is well established in the case of spectra of one-dimensional disordered 
systems: details of the distribution of random masses, random site potentials, and so on, 
only enter the periodic amplitudes of Lifshitz tails 0, [1^, [1^. 

In the present case, just as in trapping problems P, H, |10[, the quantities of most 
interest just involve the constant Fourier component of the model-specific periodic am- 
plitudes. For the total number of even-visiting 4/c-step walks, the amplitude (|6.12|) has 
been determined very accurately by comparing the data of the exact enumeration proce- 
dure to the analytical asymptotic estimate (|6.11|). 



Our analytical investigations have also led to the prediction ( |7.3| ) that all the cumu- 
lants of the maximum height M reached by an ?7,-step even- visiting walk scale as n^^^, 
with known prefactors. Usual random walks (for which M/n}/"^ has a non-trivial limit 
law) and even- visiting ones (for which M/ n^/'^ becomes more and more peaked as n gets 
large) are therefore very different in that respect. We have also investigated amplitude 
ratios, associated with higher moments of the Riccati variables, which exhibit a highly 
non-trivial dependence in the variable a. 

Finally, the asymptotic result ( |6.11|) demonstrates that the outcome ( |1.2|) of Derrida's 



original argument is very accurate in the one- dimensional it just misses a power 

of the number of steps (in one dimension we have VL = 2 and j = tt/2). The analogy with 
Lifshitz tails also suggests that the situation is less under control in higher dimension: 
going beyond the leading estimate ( |1.2| ), which involves only the volume of the Lifshitz 



sphere, indeed remains a difficult open problem |16 
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